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Abstract We report on the non-perturbative computation of the running cou- 
pling of two-flavour QCD in the Schrodinger functional scheme. The corre- 
sponding A-parameter, which describes the coupling strength at high energy, is 
related to a low energy scale which still remains to be connected to a hadronic 
"experimentally" observable quantity. We find the non-perturbative evolution 
of the coupling important to eliminate a significant contribution to the total 
error in the estimated A-parameter. 
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1 Introduction 



Under lattice regularization predictions of renormalized quantum field theories 
emerge as universal properties of critical points of models in the appropriate 
universality class. In this way the theory is defined independently of perturba- 
tion theory and may for instance be evaluated numerically. Predictive power 
resides in a surplus of relations between observables over free parameters in 
the model, and it becomes a well-defined question which part of these relations 
is amenable to approximation by renormalized perturbation theory. In QCD 
the standard expectation is that quantities associated with energies large com- 
pared to typical hadron masses can be perturbatively related to each other. 
If one is limited to this calculational framework, a small number of input 
parameters associated with large normalization energy fi, like the coupling 
a s (fi) and the quark masses for each flavour ?Sf(/i), have to be determined 
from experiment and then lead to many successful predictions of perturbative 
QCD. 

By lattice techniques it becomes possible to look beyond the perturbative 
horizon. Consequently a lot of activity goes and went into extracting infor- 
mation on the hadronic low energy sector. In particular the free parameters 
are determined in this case by inputting quantities like some hadron masses 
or Fn . Then the high energy sector can in principle be predicted by evaluating 
ct s (fJ>) and the quark masses for \i S> F n . Such calculations relating differ- 
ent orders of magnitude of physical scales represent a formidable numerical 
problem. Beside the dissimilar physical scales, infrared and ultraviolet cutoffs 
have to be extrapolated to their respective limits. A number of such direct 
approaches have nevertheless been tried, and some have found their entry into 
the particle data table as one of the most accurate determinations of 

a(Mz). In view of the very limited parameter range accessible to simula- 
tion we find it difficult to be confident about the systematic errors of these 
determinations. A computation requiring more steps but also offering much 
more control of systematic errors becomes feasible by the recursive finite size 
method using the Schrodinger functional. This technique has been developed 
by our collaboration over the last years and is reviewed in || and ||. 

In this publication we present first numerical results toward the extension 
of the method from the zero flavour (quenched) approximation to QCD with 
two light flavours which are taken massless. Sect. 2 summarizes the most 
essential results of our previous work in the present context. In sect. 3 the 
computational strategy for the A-parameter characterizing the coupling at 
large energy is outlined, followed by numerical results in sect. 4 and some 
conclusions. 
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2 Schrodinger functional setup 



To connect hadronic and perturbative scales in QCD an intermediate renor- 
malization scheme has been devised where the finite system size L is used as a 
renormalization scale. More precisely, we consider the Schrodinger functional 
given by the partition function of QCD on a cylinder of size TxL 3 in euclidean 
space 

„-r _ f TVrr „/. .11 „-S 



D[U^,^]e- b . (2.1) 

TxL 3 

In the lattice regularized form we integrate over SU(3) gauge fields U with 
the Wilson action and two flavours of 0(a) improved Wilson quarks 
Boundary conditions in the spatial directions of length L are periodic Q for 
U and periodic up to a global phase || 6 = ir/5 for tp, -ip. In Euclidean time 
Dirichlet boundary conditions are imposed at xq = 0, T by fixing spatial links 
to diagonal SU(3) matrices that are precisely specified in terms of L and two 
angles rj and v (point 'A' in and we also take T = L. The quark fields 
on the boundary surfaces jlO| are given by Grassmann values p,p and p',p', 
which are used as sources that are set to zero after differentiation. 

To achieve the convergence to the continuum limit at a rate proportional 
to the squared lattice spacing a 2 a number of irrelevant operators have to be 



tuned. The coefficient csw of the clover term [ 11 1 is set to the non-perturbative 



values quoted in parameterized form in [12]. In the Schrodinger functional at 
vanishing quark mass, that we consider here, the coefficients Ct and c*t of two 
additional boundary counter terms [1C] have to be adjusted. Here we have to 



content ourselves with perturbative estimates at one and two-loop accuracy 

MM- 



Since for Wilson fermions chiral symmetry only emerges in the continuum 
limit, the bare mass parameter is additively renormalized. For this reason we 
trade it for a quark mass defined by the PC AC relation evaluated using suitable 



states p7 [. Let /a(^o) an d fp(xo) be the matrix elements of the axial current 



and the pseudoscalar density defined in (2.1) and (2.2) of 17] with the gluonic 



boundary fields assuming the values quoted above. We form the ratio 

{ v U^o + d* )f A (x Q ) + c A ad d* fp(x Q ) 

m{x ) = . r (2.2 

2/p(sq) 

with forward (backward) derivative d (9g). For the current improvement 
coefficient ca its one-loop value [^] is taken. We now define the bare current 
mass 

J m(T/2) for even T/a , . 

mi = \ [m((T - a) /2) + m((T + o)/2)]/2 for odd T/a. ( ' 
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An alternative definition m<i just differs by averaging m(xo) and m'(xo), where 



the latter is defined [17] with the sources p', p' at the xq = T boundary leading 
t° /a'/p- These masses are expected to differ at 0(a 2 ). With either of them 
vanishing, the chirally symmetric continuum limit may be approached. 

The coupling g 2 and the additional universal dimensionless observable v 
are related to the Schrodinger functional by 

? = (2-4) 



di] \g 2 {L) 

where k is a known ||] normalization fixed by demanding g 2 = c7q + 0(<7q) with 
the bare coupling g$. 



3 Computational strategy for the A-parameter 



Our method to extract A, which characterizes the behaviour of g 2 at asymp- 



totically large energy, follows the strategy used in 18 j. By continuum extrap- 
olation we construct the non-perturbative step scaling function (SSF) 



a{u) = 9 2 (2L)\- g2{L)=umi 



(3.1) 



for a number of u- values such that by interpolation we control it over the range 
that will be needed. Then a value u max is selected (initially by guesswork) such 
that the associated scale L max where g 2 (L max ) = u max is in the hadronic range. 
By recursively solving n times 



a{g 2 {L/2))=g 2 {L) 



(3.2) 



starting with L = L m3iX we obtain values for g 2 (2 n L max ). Finally, for suffi- 
ciently large n, this coupling is perturbative and we use 



AL r 



2 n (b g 2 



x exp 



-6i/26g 



exp 



1 



dx 



2b g 2 
1 



+ 



b 2 x 



(3.3) 



_P(x) ' box 3 

to derive A in terms of L max . Here the three- loop /?- function]] for the SF- 



scheme with two flavours [14| is used and bo,b\ are its universal coefficients. 
On the right-hand side g 2 is understood to be inserted at the scale 2 _ ™L max . 

We now use 6 2 = 0.06/(4tt) 3 given in the second erratum to which has become 
necessary due to the revision of 113] in May 2001. Since our analysis depends on these 
results, an independent check seems desirable. It will be partially supplied in the near future 
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The admissibility and accuracy of renormalized perturbation theory can be 
probed by checking the stability of the result when varying n. 

In a later series of simulations we shall have to relate L max to a truly 
physical scale, for instance by computing L^^F^. The guess for u max will 
be confirmed then, if a number of order one is found, i. e. the multiple scale 
problem is avoided. The relation between A, which corresponds to the SF- 
scheme, and the MS-scheme is given by ]8|,|l4|] 

A^g = 2.382035(3)A. (3.4) 



4 Numerical results 



The continuum SSF is given by the limit 



a(u) = lim T,(u,a/L), (4.1) 

a/L^Q 



where £ is defined like a in eq. ( |3.1|) but interpreted at finite resolution a/L. 
Since PCAC becomes an operator relation in the continuum limit only, we 
adopt the convention to always tune m\{L/a) to zero on the small lattice. 
The corresponding value m\(2L/a) measured at resolution a/2L is expected 
to differ by 0(a 2 ) from nii(L/a). In the same way we expect mi(L/a) — 
TO2 (L/ a ) = 0(a 2 ) for our alternative definition m 2 of the bare current quark 
mass. We tested these expectations on our data and compared with one-loop 
perturbation theory in Fig. 1 for several couplings. Where available we include 
together with our present iVf = 2 data also quenched and bermion |l9| (iVf = 
—2) results. We conclude that lattice artefacts behave non-pathologically and 
similar to perturbative expectations. For the accessible range of resolutions 
they happen to be dominated by terms of higher order than the expected 
a 2 -contributions. 

Another place to study the approach to universal continuum behaviour is 
the relation between v and g defined in ([2.4j), 

v = uj(g 2 )= lim n(g 2 ,a/L). (4.2) 

a/L^O 

In perturbation theory Q is known to two- loop order, 

n(u, a/L) = fa + v 2 u) [l + ei(a/L) + e 2 (a/L)u + 0(u 2 ), (4.3) 

and ei,e2 encode the perturbative artefacts. In Fig. 2 data for two different 
values of the coupling are plotted. The square symbols refer to the improved 
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Figure l: Lattice artefacts of PC AC masses. Perturbative results for integer 
L/a are connected by lines. 

observable, where the Monte Carlo results have been divided by the pertur- 



bative lattice artefacts 1 + e\(a/L) + e2{a/L)u as first discussed in |2C]. We 
conclude that a much smoother continuum limit is achieved in this way. 

At present, due to limitations of computing power, we have results with 
sufficient statistics only for g 2 on lattices with L/a < 12 . The resulting values 
of £ are collected in table 1. The algorithmic aspects of these simulations have 



been presented in [21|. 



To estimate a continuum value for a from lattices with L/a = 4,5,6 
(together with the lattices at the doubled lengths) we adopt the following 
procedure. First we perturbatively correct the data with a factor analogous 
to the one in eq. ([4.3[) 



£(u, a/L) -» £( 2 > (u, a/L) = - - . /rN 2 (4.4) 

v ' ' K 1 ' 1 + 5i{a/L)u + 5 2 {a/L)u 2 v ; 

with the series for the artefacts known up to two-loop order. They depend, 
of course, on the details of the action chosen. As the two-loop boundary 
improvement coefficient Ct(go) became available only during our simulations 
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Figure 2: Circular symbols are data for v including a continuum extrapolation. 
Square symbols have perturbative lattice artefacts cancelled. 

they were partly carried out with its one- loop value (left part of table 1) and 
only later with the two-loop value. Hence two different sets of 62 had to be 
used in ( |4.4j ). We found the values of £( 2 )(u, a/L) for L/a = 5,6 constant 
within errors and fitted them to a constant (i. e. just combined them) as our 
present continuum estimates. They are found in table 2 together with £^ 2 ) 
at resolution L/a = 4 for the estimation of systematic errors (see below). In 
Fig. 3 the analogous procedure can be judged in the quenched case, where 
many more data are available. The averages of the points at L/a = 5,6 
in each series lead to the dotted lines and are to be compared with the full 
extrapolation (points at a/L = 0). 

We interpolate the values of table 2 by fitting cr(u) to a sixth order poly- 
nomial with the first three coefficients constrained to their perturbative values. 
The resulting SSF is shown in Fig. 4. It differs from the quenched SSF by an 
amount predicted well by perturbation theory for weak coupling. For values 
above about 2.5 the three-loop term contributes significantly to the /3-function 
but actually enhances the growing gap between Monte Carlo results and per- 
turbation theory. 

The fitted form for a is employed to estimate AL max in the way described 
in sect. 2 starting from u max = 3.3 and from u max = 5. Statistical errors are 
obtained by propagating the errors of the primary Monte Carlo data through 
the whole analysis, and the inclusion of another parameter in the interpolating 
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L/a 


u 


S 


u 


S 


4 


0.9793(7) 


1.0643(34) 


1.5031(12) 


1.720(5) 


5 


0.9793(6) 


1.0721(39) 


1.5033(26) 


1.737(10) 


6 


0.9793(11) 


1.0802(44) 


1.5031(30) 


1.730(12) 


4 


1.1814(5) 


1.3154(55) 


2.0142(24) 


2.481(17) 


5 


1.1807(12) 


1.3287(59) 


2.0142(44) 


2.438(19) 


6 


1.1814(15) 


1.3253(67) 


2.0146(56) 


2.508(26) 


4 


1.5031(10) 


1.731(6) 


2.4792(34) 


3.251(28) 


5 


1.5031(20) 


1.758(11) 


2.4792(73) 


3.336(50) 


6 


1.5031(25) 


1.745(12) 


2.4792(82) 


3.156(55) 


4 


1.7319(11) 


2.058(7) 


3.334(11) 


5.298(85) 


5 


1.7333(32) 


2.086(21) 


3.334(15) 


5.41(12) 


6 


1.7319(34) 


2.058(20) 


3.326(20) 


5.68(13) 



Table l: Data for the lattice step scaling function X(u, a/L). For the left part 
Ct(fl'o) was set to its one- loop value, whereas the two rightmost columns have 
been obtained with the two-loop result |L4|]. 



it 


a(u) 


£( 2 )(u,l/4) 


0.9793 


1.0768(30) 


1.0686(35) 


1.1814 


1.3277(46) 


1.3199(55) 


1.5031 


1.7489(85) 


1.7332(60) 


1.7319 


2.063(15) 


2.0562(72) 


1.5031 


1.750(8) 


1.7477(56) 


2.0142 


2.494(16) 


2.535(18) 


2.4792 


3.304(38) 


3.338(28) 


3.3340 


5.65(10) 


5.491(90) 



Table 2: Numerical results for a(u). 
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Figure 3: Quenched results for £( 2 ) to illustrate our extrapolation procedure. 
Data are mostly from the literature [9,16] apart from the two finest resolutions 
at u = 2.1 which were obtained on APEmille at Zeuthen. 




Figure 4: Step scaling function for Nf = 2 and Nf = for comparison. Dashed 
lines are perturbative results from integrations with the one-loop and, hardly- 
distinguishable, with the two- and three-loop /3-function. 
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fit for cr{u) gave only negligible changes. In this way we find the numbers in 
table 3. In the columns labelled by L/a = 4 we have replaced our continuum 





g 2 (L max ) 


= 3.3 




5 2 (-^max) — ^ 


n 


continuum 


L/a = 4 


n 


continuum L/a = 4 


5 


1.84(4) 


1.80 


6 


1.24(5) 1.19 


6 


1.86(5) 


1.78 


7 


1.26(5) 1.17 


7 


1.88(6) 


1.76 


8 


1.28(6) 1.15 



Table 3: Values estimated for — ln(AL max ) for two examples of L max . 

estimates o~(u) by Y,( 2 \u, 1/4). We regard the difference between the two 
columns as our present systematic uncertainty^] and quote at the moment 

ln(AL max ) = -1.9(2) [g\L max ) = 3.3] (4.5) 
ln(AL' max ) = -1.3(2) [f(L' max ) = 5], (4.6) 

which translates into Aj^gL max = 0.36 and Aj^gL max = 0.66 with about 



20% total errors. A corresponding number in the quenched theory [18] for 
g 2 (L max ) = 3.48 is ln(AL max ) = —1.56(8) with the full continuum extrap- 
olation and ln(AL max ) = —1.47(2) under the present procedure with only 
statistical errors indicated here. 

Finally we plot the non-perturbative evolution toward high energy for 
a(fj,) = g 2 (L)/A7r (/i = 1/L) starting from g 2 = 5 in Fig. 5. Statistical errors 
and the difference between evolving with cr(u) and £( 2 )(u,l/4) are smaller 
than the symbol size. The overall scale error implied by the uncertainty in the 



start-value in eq. (4J3), which corresponds to a rigid horizontal shift of all data 
points, is not shown here. In comparing the non-perturbative results with the 
perturbative curves we emphasize that the important point to appreciate is 
that at high energies the expected perturbative behaviour for our coupling has 
been shown to set in. On the other hand the fact that the perturbative curves 
also describe the data quite well to rather low energies refers to a property of 
our particular observable and is definitely not to be interpreted as a reflection 
of some universal property of QCD couplings. 

If, instead of evolving non-perturbatively, we had used three-loop per- 



turbation theory (eq. (3^) with n = 0) directly at the largest couplings 
g 2 = 3.3 <-> a = 0.26 or g 2 = 5 <-> a = 0.40, then we would have over- 
estimated AL max by 12% and 23% respectively. This in turn translates into 

2 A somewhat smaller value would be obtained if we took the magnitude of our per- 
turbative improvement for lattice artefacts, E*- 2 '(u, 1/5) — £(«, 1/5), as an estimate of the 
systematic error. 
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Figure 5: Evolution of a = g 2 /4tt for the Schrodinger functional coupling. 

errors of 2% and 5% for a in the range where its value is close to 0.12 corre- 
sponding to the physical value of a-^ at M%. | 

5 Conclusions 

Our results demonstrate that with the generation of parallel computers being 
installed now a computation of Aj^g including two massless flavours is becom- 
ing feasible with the ALPHA techniques. This includes - as in the quenched 
case - the possibility to probe and reduce systematic errors and, in particular, 
lattice spacing effects. Due to the high cost of the simulations, it is mandatory 
to smoothen the continuum limit as far as possible. We have therefore spent 
a significant effort on accompanying perturbative calculations. As observed 



earlier in the pure SU(2) gauge theory [20], we found that lattice artefacts 
of several quantities constructed in the Schrodinger functional are described 
quite well by perturbation theory (see Figs. 1,2). This encourages us to trust 
in perturbation theory to remove the lattice artefacts to a significant extent 
and this procedure is supported in the pure gauge theory in Fig. 3. As an esti- 
mate of the remaining systematic errors we use the difference of our results on 



3 For the estimate in the quenched approximation mentioned above the analogous error is 
smaller, since there the 3-loop /3-function happens to be closer to the non-perturbative rate 
of evolution over the relevant range. 
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the finer lattices to those on the coarsest one. Comparing with the quenched 



theory, where a robust continuum extrapolation could be carried out |jl8|, 
this error appears safe but also not over-pessimistic. Nevertheless our results 
still need to be corroborated by simulations closer to the continuum limit. We 
are in the process of simulating up to L/a = 16 to both reduce our errors for 
AL max and to put them on even firmer grounds, which will however still take 
some time. 

Already now we have clearly observed the small iVf-dependence of our 
discrete version of the j3 function (Fig. 4). For weak couplings its magnitude 
is accurately predicted by perturbation theory, while for our largest coupling 
(a ~ 0.25) it overestimates the effect significantly. In particular, the use of 
perturbation theory for couplings up to a ~ 0.4 in estimating A, would lead 
to a significant error already at our present level of accuracy. Moreover, this 
error could hardly be quantified within the framework of perturbation theory, 
which appears rather well behaved when looked at in isolation. 

Our low energy scale still has to be gauged by an experimentally observ- 
able quantity, probably by computing L^xF-ir- Also the extension to Nf = 2 
of the non-perturbative renormalization of quark masses along the lines of 
ref. [18] is within reach, once the scale dependence of a is known. In the more 
distant future we would like to include the influence of further flavours and 
their masses on the evolution of the coupling. 
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